Transverse energy fluctuations and the pattern of J/ip suppression in Pb— Pb collisions 
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The NA50 collaboration has recently observed that the J/ip production rate in Pb-Pb collisions 
decreases more rapidly as a function of the transverse energy for the most central collisions than 
for less central ones. We show that this phenomenon can be understood as an effect of transverse 
energy fluctuations in central collisions. A good fit of the data is obtained using a model which 
relates J/ip suppression to the local energy density. Our results suggest that the J/ip is completely 
suppressed at the highest densities achieved in Pb-Pb collisions. 



Among the various particles which are produced in 
nucleus-nucleus collisions at high energy, the J/tp me- 
son plays a special role. Because of the large mass of 
the charm quark, cc pairs are produced on a short time 
scale and their evolution probes the state of matter in the 
early stages of the collisions. In particular, their binding 
into J/ip mesons may be hindered by the presence of a 
quark-gluon plasma in the collision zone where they are 
produced, leading to the so-called tl J/tp suppression" |Q. 
Indeed, the rate of J/ip production in nucleus-nucleus 
collisions is smaller than what can be expected on the 
basis of extrapolations based on independent nucleon- 
nucleon collisions. However, the effect observed in col- 
lisions involving oxygen and sulphur projectiles can be 
attributed to the so-called nuclear absorption, a mecha- 
nism also at work in proton-nucleus collisions ■ But 
nuclear absorption alone is insufficient to account for the 
large suppression observed in Pb-Pb collisions, which has 
been qualified for this reason as "anomalous suppression" 
||. Further final state interactions seem to be needed 
to explain the data, and various scenarii based either 
on hadronic interactions ||, or quark-gluon plasma for- 
mation [7j|] have been proposed with various degrees of 
success 

The latest data obtained by the NA50 collaboration 
p0| allow further progress. These data provide new in- 
formation on the dependence of the J / ip production rate 
on the total transverse energy (Et). In particular, it is 
found that the rate of J/ip production decreases faster 
with increasing Et in the most central Pb-Pb collisions 
than in the less central ones. It has been suggested that 
the resulting pattern of J/ip suppression as a function 
of Et exhibits the two drops corresponding to the suc- 
cessive meltings of charmonium bound states, first the 
X (which is expected to produce 40% of observed J/ip) 
and then the J/ip ■ We shall in fact argue that the 

extra suppression observed at high Et may have a sim- 
pler origin: the largest transverse energies are achieved 
through fluctuations rather than by a change of the col- 
lision geometry, and, as anticipated in Refs. [plJl^], this 
is enough to produce an extra suppression with the right 
order of magnitude. 



The model developed in Ref. g , on which the present 
analysis relies, relates the anomalous J/ty suppression to 
the local energy density: if the energy density at the point 
where the J/ip is produced exceeds a critical value e c , the 
J/ip disappears. This model was motivated by the simple 
observation that the local energy density is higher in Pb- 
Pb collisions than in any other system which had been 
studied previously (e.g. S-U collisions), even though the 
average density in Pb-Pb collisions does not exceed that 
in other systems. This simple picture was shown to ac- 
count quantitatively for the data with a single parame- 
ter, e c . Here, we extend this model taking into account 
transverse energy fluctuations, which were neglected in 
the numerical estimates presented in || . 

The transverse energy produced in a nucleus- nucleus 
collision can be used as a measure of the impact param- 
eter, the largest values of Et corresponding to the most 
violent collisions at small impact parameter. However 
the correlation between b and Et is not one-to-one: for a 
given impact parameter, the produced transverse energy 
fluctuates. The effect of these fluctuations is particularly 
visible for the collisions producing the largest transverse 
energies: these correspond essentially to central collisions 
with nearly vanishing impact parameters. If one assumes 
that a fluctuation in Et results from a fluctuation in the 
energy density over the entire collision zone, one sees that 
an Et fluctuation is accompanied by an increase of the 
region where the energy density exceeds e c and hence can 
be responsible for an amplification of the J/ip suppres- 
sion. 

As a first estimate of this effect, let us consider central 
collisions of two identical nuclei described by sharp sphere 
densities, and ignore the nuclear absorption. At zero 
impact parameter, the energy density at point s in the 
transverse plane (see Fig.|l|) is proportional to the number 
of participants in a small tube centered at s (see below) , 
that is e(s) oc \f R 2 — s 2 . We assume that all the J/ip's 
emerging at points s < s c where the energy density is 
e(s) > e c are suppressed. The radius s c of the total 
suppression zone is then related to e c by: 
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where e max — e(s = 0). The density of cc pairs pro- 
duced at point s is proportional to the number of nucleon- 
nucleon collisions in a small tube centered at s, i.e. to 
R 2 — s 2 . The survival probability follows from a simple 



geometrical counting: 
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This result [ ft2[ provides a simple geometrical interpreta- 
tion of the effect of a transverse energy fluctuation: as- 
suming such a fluctuation to be evenly distributed over 
the whole transverse plane, an increase of Et results in 
an increase of the volume in which the J/^'s are sup- 
pressed. At the same time, the number of cc pairs, pro- 
portional to the number of nucleon-nucleon collisions, re- 
mains constant, independent of Et- This is the origin, 
in this picture, of the expected extra suppression at large 
Et- A fluctuation of Et by about 10% increases e max by 
the same amount, and, according to Eq. (j2j) , decreases the 
J/ip production by about 30%; this is indeed the order 
of magnitude of the observed effect. 




FIG. 1. Geometry of a central collision: the shaded area 
denotes the region where anomalous suppression takes place. 

We now turn to a more quantitative calculation. We 
assume that the transverse energy produced in a nucleus- 
nucleus (A-B) collision is proportional to the number of 
participants at impact parameter 6, N p (b), as predicted 
in the wounded nucleon model Jl3| . The actual depen- 
dence may be somewhat stronger than a simple linear 
relation |jlj,^5| , and such deviations are important in co- 
mover scenarios 0, but we do not expect them to alter 
our results significantly. Furthermore, as in || we make 
this relation local, and take the energy density at point 
s and impact parameter b to be proportional to the den- 
sity n p (s,b) of participants in a plane orthogonal to the 
collision axis: 



n p (s, b) = T A (s) 
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where ctjv = 32 mb is the nucleon-nucleon inelastic 
cross section, and T A (s) = J_°° dz pa(s,z). The to- 
tal number of participants at impact parameter b is 
N p (b) = J d 2 s n p (s,b). For our numerical estimates, 
we parametrize p A (r) by: 



PA (r) = 



Po 



l + exp(^) 



J p A (r)d 3 r = A, 



(4) 



with a = 0.53 fm, R A = 1.1 A 1 / 3 = 6.52 fm and po = 0.17 
fm" 3 for 208 Pb. 

At fixed impact parameter, the number of participants 
may fluctuate. We shall ignore these fluctuations, fo- 
cussing here only on the fluctuations of the transverse 
energy produced by a fixed number of participants or, 
equivalently, at a fixed impact parameter. The corre- 
sponding distribution is chosen to be the following gaus- 
sian distribution |Bm: 
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The mean value of the transverse energy is (i/r)(b) = 
qN p (b), with q the average transverse energy per par- 
ticipant, and the dispersion is <j\ = aq 2 N p (b), with 
a a dimensionless parameter. The fit to the minimum 
bias £r-distribution by NA50 yields q = 0.274 GeV and 
a = 1.27 0. Then the "knee" of the ^-distribution, 
defined as E T nee = qN p (0), sits at about 108 GeV. 

We assume, as in the simple model discussed above, 
that the fluctuation in Et at given b is distributed over 
the entire collision zone proportionally to n p (s, b) so that 
the energy density can be written as: 
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where (ET)(b) oc N p (b) is the average transverse energy. 

The production of cc pairs in an A-B collision goes like 
the number of nucleon-nucleon collisions, and so does 
that of Drell-Yan pairs. Thus the probability P(DY\b) 
that a Drell-Yan pair be produced given that a nucleus- 
nucleus collision has taken place at impact parameter b 
is given by: 
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1 — e ~ CTjvTA£! ( b ) is the probability to have at least one in- 
elastic collision (minimum bias) at impact parameter b. 
For large impact parameters, P(DY\b) becomes indepen- 
dent of b, P(DY\b) — > l G Ni which is the probabil- 
ity to produce a Drell-Yan pair in an inelastic nucleon- 
nucleon collision. 

Note that for a fixed impact parameter, the number of 
nucleon-nucleon collisions may also fluctuate. We shall 
ignore such fluctuations, recognizing that this is a source 
of uncertainty on the Et distribution of the Drell-Yan 
pair production. With this simplification, the number 
of Drell-Yan pairs is entirely determined by the nuclear 
geometry, i.e., by the impact parameter of the collision. 
Besides, the production of a Drcll Yan pair is not af- 
fected by final state interactions, thus it is not directly 
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sensitive to the transverse energy produced in the colli- 
sion (it will become indirectly related to Et through the 
relation between b and Et)- 

Final state interactions, however, strongly modify the 
J I tp production cross section. In addition to the standard 
nuclear absorption, we model the anomalous suppression 
by simply assuming that any J/ip emerging at point 
s is suppressed if the energy density e(s) is higher than 
e c . Then, the probability to produce a J /tp given that a 
collision has taken place at impact parameter b and has 
produced a transverse energy Et is written as follows: 
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where n c is a parameter proportional to the critical en- 
ergy density e c defined above. In this formula, a a ~ 6.4 
mb is the absorption cross section p0| , and a^ N is the 
J I ip production cross section in a nucleon-nucleon colli- 
sion. 

Putting everything together, we can write the proba- 
bility to produce a J jip in a collision at a given Et as: 
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By replacing in this formula P(ijj\ET,b) by P(DY\b), 
we obtain the probability to produce a Drell-Yan pair at 
a given Et- We can then construct the ratio: 
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which we shall use in order to compare the model with 
experimental data. 

Indeed, the NA50 collaboration presents its re- 
sults essentially as the same ratio, where the numerator 
is an experimentally measured quantity, while the de- 
nominator, which is model dependent, is deduced from a 
theoretical analysis that is identical to the one presented 
above. In order to complete our theoretical estimate, we 
need an extra input, the normalisation factor a^ N / a^y 
(multiplied by the branching ratio to the dimuon decay 
channel) which is obtained from experiments with lighter 
projectiles and targets Jl9| . The value of this normaliza- 
tion factor given by NA50 is a$ N /a%$ = 53.5 ± 3 @- 
The fits in Figs.|] and || have been obtained with the val- 
ues 52 and 53.5 respectively. 

A plot of the calculated AT(Et) is shown in Fig.|| to- 
gether with data from NA50. A reasonable fit of the 1998 
data is obtained with a single parameter n c = 3.7 fm~ 2 . 
Above 40 GeV, the curve deviates from nuclear absorp- 
tion because the anomalous suppression mechanism sets 
in. Note that the points at low Et from 1996 data p8| , 
on the other hand, are above the curve, i.e. they show 



less J/ip absorption than expected from an extrapola- 
tion of proton-nucleus and nucleus-nucleus data, a fact 
still unexplained. Above 100 GeV, i.e. approximately at 
the position of the knee, a second drop occurs which, as 
explained above, is associated with the increase of Et 
due to fluctuations. The fact that we reproduce quan- 
titatively the relative magnitude of this drop suggests, 
according to the discussion following Eq.(2), that all the 
J/ip mesons disappear at the highest densities achieved 
in the system. 
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FIG. 2. The J/ip survival probability in a Pb-Pb collision 
as a function of the transverse energy in GeV, measured by 
NA50 in 1996 (open circles) and in 1998 (close circles); after 
absorption in nuclear matter alone (dot-dashes); and after 
dissolution in a quark-gluon plasma with a critical density 
n c — 3.7 fm -2 , without (dashes) and with fluctuations (full 
curve) . 



We have to point out that the rapidity window used 
to mesure the J/ip and the one used to determine the 
-E/r-distribution are different. Indeed, muons pairs are 
selected for 2.82 < y < 3.92, whereas the neutral trans- 
verse energy is measured in the 1.1 — 2.3 pseudo-rapidity 
window |18|. In the calculation presented above, we have 
implicitly assumed that the fluctuations in these two ra- 
pidity windows are strongly correlated. 

If, on the contrary, there were no correlation between 
Et fluctuations and the energy density at the point where 
the cc pair is produced, then the suppression criterion 
should involve only the average transverse energy at a 
given impact parameter, and (||) should be replaced by 
n p (s, b) > n c . The resulting M(Et) is shown as the 
dashed curve in Fig.|. In this scenario, the suppression 
is a function of the impact parameter b only: as a con- 
sequence, the suppression saturates at large Et, to its 
value at b = 0. Such a model for J/ty suppression is pro- 
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posed in 0,0, where, indeed, the criterion for suppres- 
sion depends only on b. Although it refers to the density 
of strings rather than to that of participants, this bears 
little influence on the results. 

Clearly, our theoretical curve in Fig.|^ does not provide 
a perfect fit to NA50 data, which have small error bars. 
In order to improve the quality of the fit, we did an- 
other calculation with a more gradual suppression mech- 
anism: namely, we replaced the 0-function in Eq.Q by 
[1 — tanhA(n — n c )]/2. The results for n c = 3.75 fm -2 
and A = 2 fm 2 are displayed in Fig.[| A perfect fit of the 
1998 data is obtained over the entire Et range from 40 
GeV to 120 GeV. Note that above 80 GeV, the suppres- 
sion becomes total at the hottest point of the system, so 
that the behaviour at large Et is essentially the same 
as in Fig.||. We obtain a fit indistinguishable from that 
of Fig.|| with a two-threshold scenario |l(J, by assum- 
ing that 40% of the J ftp are suppressed if the density 
lies between n c i — 3.3 fm~ 2 and n c2 = 4.0 fm~ 2 and 
100% above n c2 . Note that the second threshold does 
not produce any visible structure in the Et dependence 
of a^/any 



In summary, we have shown that the second drop in 
the J ftp yield around the knee of the £/r-distribution 
can be interpreted as an effect of Et fluctuations, if one 
assumes that J ftp suppression is related to the local en- 
ergy density in the system. We have achieved our best 
fit by assuming that the suppression increases gradually 
with the energy density; in particular, there is no indi- 
cation in the data that the suppression occurs in two 
steps. However, in order to reproduce the magnitude of 
the observed "second drop" , we need to assume that all 
the J ftp's disappear at the highest densities achieved in 
a Pb-Pb collision. 
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FIG. 3. Same as previous figure, with a gradual onset of 
anomalous J ftp suppression (see text). 

Finally, we would like to compare our model with that 
recently proposed in ^0|, where the J/ty suppression is 
explained by interactions with comoving particles. The 
density of comovers is related to the transverse energy 
in the same way as the energy density in our model. 
Thus these authors also obtain increased J/^> suppres- 
sion, but much less pronounced than in NA50 data. This 
is because, as explained in detail in jl2|, absorption by 
comovers depends more weakly on tranverse energy fluc- 
tuations than in the present model. 
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